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Abstract. We analyze the contribution of gamma-ray pulsars from the first Fermi- 
Large Area Telescope (LAT) catalogue to the local flux of cosmic-ray electrons and 
positrons (e + e _ ). We present new distance estimates for all Fermi gamma-ray pulsars, 
based on the measured gamma-ray flux and pulse shape. We then estimate the 
contribution of gamma-ray pulsars to the local e + e~ flux, in the context of a simple 
model for the pulsar e + e~ emission. We find that 10 of the Fermi pulsars potentially 
contribute significantly to the measured e + e~ flux in the energy range between 100 
GeV and 1 TeV. Of the 10 pulsars, 2 are old EGRET gamma-ray pulsars, 2 pulsars 
were discovered with radio ephemerides, and 6 were discovered with the Fermi pulsar 
blind-search campaign. We argue that known radio pulsars fall in regions of parameter 
space where the e + e" contribution is predicted to be typically much smaller than from 
those regions where Fermi-LAT pulsars exist. However, comparing the Fermi gamma- 
ray flux sensitivity to the regions of pulsar parameter space where a significant e + e~ 
contribution is predicted, we find that a few known radio pulsars that have not yet 
been detected by Fermi can also significantly contribute to the local e + e~ flux if (i) 
they are closer than 2 kpc, and if (ii) they have a characteristic age on the order of 
one mega-year. 
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1. Introduction 

In the last two years, data on high-energy cosmic-ray electrons and positrons have 
attracted a great deal of interest. The positron fraction measurement reported by the 
PAMELA Collaboration in Ref. [1], combined with the high-statistics electron-positron 
(e + e~ ) flux measurement by the Fermi Large Area Telescope (LAT) Collaboration [2] 
are widely believed to conclusively point to a new class of primary positron sources in the 
Galaxy (see e.g. [3]). Interestingly, one such class that has been widely advocated and 
discussed of late is the pair-annihilation or the decay of Galactic particle dark matter 
(see e.g. [4]). A more mundane possibility, but by no means more or less disfavored 
by currently available data, is that the new class of positron sources is associated to 
nearby mature rotation-powered pulsars, as earlier envisioned in relation to the recent 
cosmic-ray data e.g. in Ref. [5-9]. 

While from the point of view of advancing our understanding of the fundamental 
structure of particle physics the dark matter scenario appears more appealing, pulsars 
might constitute if not the entirety of the detected positron excess, at least a potentially 
very significant background to searches for new physics. In this respect, gaining an 
accurate understanding of the contribution of pulsars to the local flux of high-energy 
e + e~ addresses questions that go beyond the realm of classical high-energy astrophysics 
and that are of relevance to the quest for physics beyond the Standard Model as well. 

The Fermi-LAT Collaboration recently delivered [10] the first catalogue of 46 high- 
confidence gamma-ray pulsars, based on the first six months of LAT data. Roughly half 
of these pulsars were discovered using radio ephemerides, i.e. monitoring the locations 
and timing, in the gamma-ray sky, of known radio pulsars. The remaining pulsars 
were discovered by looking for a pulsed emission in bright gamma-ray sources, either 
corresponding to the potential locations of neutron stars identified via observations at 
other wavelengths, or via the pulsar blind-search campaign (see Ref. [11-13]); finally, 6 
pulsars correspond to previously known EGRET gamma-ray pulsars [14, 15]. 

The scope of the present study is to assess the contribution of pulsars in the 
first Fermi-LAT gamma-ray pulsar catalogue to the local flux of high-energy cosmic- 
ray electrons and positrons. There are several reasons to believe that gamma-ray 
pulsars contribute significantly to the local e + e~ flux. On theoretical grounds, the 
mechanism through which gamma rays are produced in the pulsar magnetosphere 
necessarily requires the generation of high-energy electron-positron pairs, independent 
of whether particle acceleration and radiation occurs in the neutron star polar cap (see 
e.g. Ref. [16-21]) or in the so-called outer gap ([22-24]). A possible caveat is the case of 
gamma rays produced by curvature emission from primary out-flowing electrons escaping 
magnetized neutron stars: in this case no positron counterpart would be produced by 
the gamma-ray emitting object. 

While over 1,800 radio pulsars are known [25], it is widely believed that these 
radio-loud objects only constitute a fraction of the wider class of rotation powered 
pulsars. The pulsar radio catalogue completeness is not only flux-limited, but - more 
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importantly - it is limited by the pulsar radio beam geometry: a radio pulsar is visible 
only if the relatively narrow radio beam crosses the observer's line of sight. Gamma-ray 
pulsars, on the other hand, are believed to exhibit much wider beams: powerful nearby 
e + e~ injectors might thus very well be radio-quiet due to geometry, and yet produce a 
significant flux of gamma radiation. An outstanding example is the Geminga pulsar, 
that has been advocated e.g. in Ref. [5, 6] as the possible main contributor to the excess 
of cosmic-ray positrons detected by PAMELA. We therefore find it compelling to use 
the recent high-confidence list of gamma-ray pulsars to study the contribution of this 
class of primary positron emitters to the local e + e _ flux. 

One of the main difficulties in assessing the contribution of gamma-ray pulsars to 
the local e + e~ cosmic-ray flux is the lack of information on the pulsar distance, especially 
if the pulsar has not been detected yet at other wavelengths. In the present analysis, 
we start, in Sec. 2, from a systematic evaluation of the distance to the relevant gamma- 
ray pulsars in the first Fermi-LAT pulsar catalogue, calibrated on previously available 
information on distances to selected objects, and making use of information from the 
measured gamma-ray flux and pulse shape and from the simulation results of the gamma- 
ray pulsar "Atlas" of Ref. [26] . We then proceed in Sec. 3 to evaluate the contribution 
from the Fermi-LAT gamma-ray pulsars to the flux of local leptonic cosmic rays, in the 
context of a simplified model for the flux of e + e~ injected in the ISM from the pulsar 
magnetosphere, and we compare with existing recent related studies [5-9]. We then 
explore the relevant pulsar parameter space, with respect to the expected contribution 
from radio-quiet versus radio- loud pulsars (Sec. 4), summarize and conclude (Sec. 5). 

2. Gamma-ray pulsars: Distance Estimates 

A pulsar is believed to consist of a rapidly rotating neutron star radiating across different 
wavelengths. Pulsars are remarkably stable clocks, emitting radiation with very precise 
time periods, with little or no variablity. While first discovered at radio wavelengths, 
pulsars have been detected across all wavelengths, including in the high-energy regime. 
With the launch of the space-based gamma-ray telescopes EGRET and of the current 
generation instrument, the LAT on board the Fermi Space Telescope, pulsars have been 
discovered to pulse in gamma-rays, a crucial step to pinpoint the mechanism itself by 
which pulsars radiate. 

Gamma-ray pulsars are rotation-powered pulsars, where the electromagnetic 
radiation is fueled through the loss of angular momentum. While only a tiny (of the 
order of 10~ 6 ) fraction of rotational energy is emitted at radio frequencies, the fraction 
of spin-down luminosity converted into gamma-rays is typically much larger, and not 
infrequently of order unity, making pulsar rotation directly responsible for their high- 
energy radiation (see e.g. Ref. [10,26]). Broadly speaking, gamma-ray pulsars exist 
in two distinct populations: young, energetic pulsars, and old, recycled millisecond 
pulsars. Since e + e~ are injected in the ISM in the early stage of a pulsar's lifetime, 
and since the very efficient energy losses for high-energy e + e~ limits their injection age, 



Gendelev, Profumo & Dormody: e + e from Fermi Gamma-Ray Pulsars 



4 



millisecond pulsars are not relevant for our present purposes, and we thus only focus 
here on energetic, relatively young (characteristic age less than 10 Myr) pulsars. These 
pulsars spin with typical frequencies in the range 0.25-100 Hz and feature relatively large 
spin-down energies E = 10 34 — 10 36 erg s -1 . Electrons and positrons produced in the 
pulsars' magnetosphere, however, are trapped by the termination shocks of the pulsar 
wind neabula and flow into the ISM only when the envelope is suppressed enough. This 
process of diffusion of the e + e~ produced in the regions around the neutron star into 
the ISM has a characteristic timescale ranging between 10 4 and 10 5 yr (see e.g. [27]). 
Cosmic-ray e + e~ from very young pulsars, like Vela, might therefore not have diffused 
into the ISM yet [9] and should not contribute significantly to the detected e + e~ local 
flux. The timescale of e + e~ injection in the ISM therefore generically further cuts out 
the youngest pulsars as contributors to the local e + e~ flux. 

The detection of gamma-ray pulsars proceeds via two possible pathways [10]: (i) by 
using the localization and timing information provided by radio ephemerides and folding 
on the gamma-ray photon data to detect gamma-ray pulsations (known as radio-selected 
gamma-ray pulsars), and (ii) by calculating periodicity in the gamma-ray photon data 
alone using blind frequency searches (known as gamma- selected gamma-ray pulsars [12]). 
Both methods are robust and have yielded the discovery of close to 50 high-significance 
gamma-ray pulsars in total [10]. 

The spin-down luminosity E of a rotation powered pulsar is simply the derivative 
of the rotational kinetic energy E = \Ioj 2 , E = Iuu. The characteristic age t c h is the 
time it would take an ideal magnetic dipole born spinning at infinite frequency to slow 
to the observed pulsar frequency today. It is only a rough estimate of the actual age 
of the pulsar, but it provides useful information in determining which population of 
pulsars might contribute the most to the e + e~ local flux, since as we shall see below 
the injection time crucially enters the diffusive propagation of charged leptonic cosmic 
rays. The characteristic age can be calculated simply from the spin parameters uj and 
w: t c h = —w/ (2w). 

To calculate the distances of the gamma-ray pulsars detected by Fermi-LAT, we 
use the heuristic model of Ref. [26], 

( E \ V2 

L^wE^Cx — x 10 33 ergs' 1 , (1) 

\ 10^ erg s~ x / 

where L 1 is the gamma-ray luminosity in the 0.1 GeV to 100 GeV energy range, E is the 
spin-down luminosity of the pulsar, w oc E~ x l 2 is the gamma-ray efficiency, and C is a 
constant of order unity that is in principle theoretically calculable given a pulsar gamma- 
ray production mechanism. Here, we take a phenomenological approach and we calibrate 
the coefficient C by fitting it to Fermi-LAT gamma-ray data for the 15 pulsars from the 
brightest LAT source list catalogue [28] with solid associated distance determinations 
(the latter taken from Ref. [10]). We also considered the theoretically naive possibility 
that the gamma-ray luminosity is simply proportional to the available rotational energy 
of the neutron star, i.e. L 1 oc E, although somewhat disfavored by data from the 
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above mentioned gamma-ray pulsars with distance determinations [10]. In this case, 
the gamma-ray luminosity of pulsars with a low spin-down luminosity is suppressed, 
leading in turn to a suppressed estimate for the distance to most of the relevant 
gamma-ray pulsars. This implies therefore that the assumption that L 7 oc E leads 
to an over-estimate of the e + e~ contribution from gamma-ray pulsars. In other words, 
the assumption we employ for the gamma-ray pulsar luminosity leads to conservative 
predictions for the flux of cosmic-ray e + e~ . 

The most precise distances to pulsars are determined from parallax measurements, 
but this method is only reliable for small distances (< 500 pc). For pulsars with no or 
unreliable parallax data, one must make use of other means of distance determination, 
such as radio Disperson Measure or X-ray absorption. In particular, the Dispersion 
Measure (DM) distance estimate hinges upon the effect on the pulsar radio beam of 
the integrated column density of free electrons between the Earth and the pulsar. The 
interaction between radio pulses and the charged particles (i.e. free electrons) in the 
Galaxy causes the pulses to be delayed by an amount depending on the electron density 
and on the radio frequency. By measuring the frequency shift and using a model for the 
distribution of free electrons, one can assess the distance to radio pulsars. Most radio 
pulsar distances in the Fermi Pulsar Catalog [10] are obtained in this way. Dispersion 
measurements in Ref. [10] make use of the NE2001 electron density distribution model 
of Ref. [29] . Other distance measurements used in the first LAT pulsar catalogue include 
using the Doppler shift of absorption or emission lines in the neutral hydrogen spectrum 
(HI), and X-ray flux measurements (see Ref. [10]). 

To estimate the distance to a pulsar with gamma-ray measurements, we use the 
relation [26] 

L 7 = f n (4nD 2 )G„ (2) 

where fn is the flux correction factor, which is a pulsar-dependent and model-dependent 
quantity that translates between the phase-averaged flux observed along the Earth line 
of sight and the gamma-ray flux averaged over the entire sky (see Ref. [30]), D is the 
estimated distance to the pulsar, and G 7 is the 7-ray energy flux inferred by the Fermi 
analysis of the gamma-ray flux and spectra [10] (this quantity is indicated in [10] with 
the symbol Gioo)- 

The coefficient C is assumed here to be universal, and to estimate its value we 
restrict our analysis to pulsars closer than 10 kpc and with a characteristic age t c h < 10 7 
yr. Incidentally, this implies that the spin-down luminosity of the pulsars we consider 
have spin-down luminosities E > 0.5 x 10 34 erg s -1 . While for very distant pulsars the 
uncertainty on the distance from radio measurements is very significant, the constraint 
on age stems from the fact that the maximal energy e + e~ produced 10 7 yr ago can have 
today is smaller than roughly 20 GeV, and therefore pulsars older than 10 7 year cannot 
contribute significantly to the local e + e~ flux in the energy range of interest here. Notice 
also that this cut on the pulsar age excludes milli-second pulsars as well. Interestingly, 
the range of spin-down luminosities selected by the cuts we employ corresponds to one 
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where the heuristic formula of Eq. (1) should hold: in several motivated pulsar gamma- 
ray emission models, including the slot and outer-gap frameworks, below a spin-down 
luminosity of 10 34 erg s -1 the gap saturates and the potential drop is almost entirely 
employed to fuel the pair cascade processes, breaking the relation in Eq. (1), as pointed 
out e.g. in Ref. [31,32]. The constraints on age and distance leave 10 pulsars out 
of the 15 LAT bright pulsars with solid distance determinations, namely J0534+2200 
(Crab), J0835-4510 (Vela), J0633+1746 (Geminga), J1709-4429, J1952+3252, J1057- 
5226, J1048-5832, J0631+1036, J1028-5819 and J1509-5850. 

The value for the flux correction factor fa is calculated for each individual pulsar, 
and derived using the results of the "Atlas" presented in [26], and making use of 
predictions from both the Outer Gap (OG) and the Two-Pole Caustic (TPC) models. 
The correction factor also depends on the magnetic inclination angle ( and on the 
observer viewing angles a; Since we do not have exact values for either one of these 
angles, we can at best find a range of values for fa from the peak separation observed and 
reported in the first LAT pulsar catalogue, Ref. [10]. The first step is to choose the most 
plausible model for fa: Following the results of Ref. [26], for those pulsars which have 
a peak separation of 0.15-0.30 the OG model is the clear choice, while for separations 
larger than 0.47 the TPC model is typically more adequate. For intermediate peak 
separations, we find that the OG model predicts values of fa giving distance estimates 
closer to the catalogue distances in almost all cases. To estimate the coefficient C we 
therefore employ, for intermediate peak separations, the fa predictions from the OG 
model. 

Given the assumptions above for the pulsar model, we then estimate fa for each of 
the 10 pulsars we use to calibrate the parameter C in the following way: we derive the 
value of the gamma-ray efficiency w from Eq. (1) and use the results of the simulations of 
[26] corresponding to the value of w closest to what we derive from the actual measured 
gamma-ray and spin-down luminosities. We then utilize the measured gamma-ray peak 
separation between the two most prominent peaks, as found in [10], to record the range 
of viewing and magnetic inclination angles. Finally, this range is used to infer, again 
using the results of the simulated pulsars of Ref. [26], the flux correction factor for 
each individual pulsar. Our final estimate for the parameter C is C = 6.04 ±0.75, where 
the quoted error results from the uncertainty on fa and on G 7 for the 10 pulsars we use 
to estimate the parameter. 

With the determination of the best fit value for C, with the method to extract fa 
outlined above, and with data on the measured gamma-ray energy fluxes and pulse 
shapes, we can obtain distance estimates for the Fermi-LAT gamma-ray pulsars of 
Ref. [10]. In turn, these distances crucially enter the estimates of the local e + e~ flux 
produced by known gamma-ray pulsars. The uncertainties associated to the distance 
estimates stem from (i) the uncertainty on the Fermi-LAT measured gamma-ray energy 
flux Gj, from (ii) the estimate of fn, the flux correction factor, and from (iii) the 
parameter C. We use standard error-propagation techniques to calculate the propagated 
uncertainty on the distance. The other quantity entering the gamma-ray pulsar distance 
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Figure 1. A comparison of the pulsar distances calculated with Fermi-LAT gamma- 
ray flux data and of distances obtained, for some of the same pulsars, with other 
methods (indicated with a blue double arrow; see Ref. [10]). Distances, in kpc, are 
plotted against the pulsars' spin-down luminosity E in erg s~ l . Black circles indicate 
pulsars discovered in the Fermi-LAT pulsar blind search [12], while red squares refer 
to either EGRET pulsars [33] or to pulsars found with radio ephemerides. 



estimates, namely the pulsar spin-down luminosity, is inferred from pulsar timing 
information and has virtually no associated uncertainty. 

Tab. 1 quotes the distances for the pulsars in the Fermi-LAT catalogue and 
compares them with our estimates based on the measured gamma-ray flux and gamma- 
ray pulse properties (the pulsar names are listed in the first column, the catalogue 
distances, when available, in the second column and our estimates in the third column). 
We include only pulsars with a characteristic age below 10 Myr (older pulsars do not 
contribute to the cosmic-ray e + e~ flux in the energy range we are interested in here) 
and nearby (i.e. distances less than 10 kpc) objects. The first ten pulsars in the table 
correspond to the top 10 pulsars contributing to the local flux of cosmic-ray e + e~ (we 
present details on how we estimate this flux in the next section). 

We show in Fig. 1 a visual comparison of the DM distances and gamma-ray 
distances, as a function of the pulsar's spin down luminosity E, for those same 10 e + e~ 
-brightest pulsars. Objects indicated in black correspond to pulsars found with gamma- 
ray data, while those in red to either previously known gamma-ray pulsars (EGRET 
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pulsars) or pulsars discovered with radio ephemerides. The thick blue double-arrowed 
lines indicate alternate distance estimates for the pulsars under considerations, for those 
pulsars with a distance determination quoted in Ref. [10]. 

We find that several of the distance estimates perfectly agree with the known DM 
pulsar distances from the catalog (e.g. J1836+5925, J1057-5226, J0633+1746, and 
J1809-2332; the distance estimates are compatible at the 2a level for J1741-2054 and 
for J2043+2740). Only in two cases do we find that the gamma-ray distance either 
overestimates (for J0659+1414) or underestimates (for J2021+4026) the DM distance 
as determined in the Fermi-LAT pulsar catalogue from the NE2001 model. In the case 
of PSR J2021+4026, the gamma-ray luminosity reported in the pulsar LAT catalogue 
[10] includes a bright un-pulsed emission, potentially associated to a secondary nearby 
source of gamma rays, or to magnetospheric emission [10]. The LAT catalogue [10] 
also points out PSR J0659+1414 as an outlier with an inferred luminosity a factor 30 
below the general trend for other gamma-ray pulsars. One possibility for the apparent 
very low gamma-ray efficiency of PSR J0659+1414 is the smallness of the viewing angle 
C < 20°, as envisioned e.g. in Ref. [34]. 

Under- and over-estimates of pulsar distances determined with gamma-ray data 
are therefore dominantly due to the pulsar-dependent gamma-ray conversion efficiency, 
which in our setup is assumed to essentially be universal and mirrored by the average 
value of the constant C (although some pulsar-dependent information on the flux 
correction factor fa is included in the present analysis). Nevertheless, the agreement 
we find between DM distances and gamma-ray distances for pulsars contributing 
significantly to the local e + e~ flux is remarkable and gives us confidence on the distance 
extrapolation for pulsars for which no DM estimate is known. Finally, we also point 
out that among pulsars that do not contribute significantly to the local e + e~ flux (lower 
part of Tab. 1), we also find very good agreement between our gamma-ray distance 
predictions and distance estimates as listed in Ref. [10]. The distance estimates we 
obtain with gamma-ray data all agree with other distance determinations at the 2-a 
level or better. 

3. The Contribution of Gamma-ray pulsars to the Cosmic-Ray e + e Flux 

The e + e~ output from pulsars has been investigated recently in several studies (e.g. 
Ref. [5-7,9]), given the interest in primary sources of energetic cosmic-ray positrons 
sparked by the PAMELA results on the positron fraction [1]. The level of emission 
of high energy e + e~ is expected to be pulsar-dependent, and no consensus exists on a 
robust theoretical prediction for the spectrum and normalization of this flux. For the 
sake of illustration, we adopt here a simple and easily reproducible model for the e + e~ 
emission from pulsars, and for the propagation of high-energy e + e~ from the neutron 
star to the Earth. Also for simplicity and to allow for a comparison among the various 
pulsars we adopt the same parameters and assumptions for all objects. We closely follow 
here the notation and the analytical expressions of Ref. [3] and references therein. 



Gendelev, Profumo & Dormody: e + e from Fermi Gamma-Ray Pulsars 



9 



We assume an e + e spectrum at source of the form 



E 



Q(E,t,r) = Q [YGeVj e M-E/E cut ]5(t - t )5(f), (3) 

and we take a spectral index T = 1.7. This value is in the range of the gamma-ray 
spectral indexes reported in the Fermi-LAT pulsar catalogue [10], where 1 < T 7 < 2 
with peaks at 7 ~ 1.3 and 1.7. Although the physical mechanism that produces e + e~ 
in the pulsar magnetosphere is presumably the same that then gives rise to the gamma- 
ray emission [35], reacceleration in the pulsar wind nebula (PWN) or in the supernova 
envelope can significantly affect the e + e~ spectral index. The value of T ~ 1.7 we 
adopt here is however consistent with the e + e _ population plausibly responsible for the 
synchrotron and inverse Compton emission from e.g. the Crab PWN [36] (see also the 
discussion in [9]). 

We adopt the value E cut = 1 TeV, compatible with the expected decrease with 
time (see e.g. Ref. [37]) of the large cutoff energy expected from the observed very 
high-energy gamma-ray emission from the PWN of young pulsars [38]. We neglect here 
the time delay between the pulsar birth and the point in time when e + e~ are injected in 
the inter-stellar medium (to in Eq. (3)) because (i) the uncertainty in the characteristic 
age versus the actual pulsar age is of the same order of magnitude as the ISM injection 
time scale, and (ii) our results are basically insensitive to this parameter for the pulsars 
that contribute the most to the local e + e" flux. 

Finally, the normalization parameter Q Q was set for each pulsar to the value such 

that 

/ Ex Q(E)dE = E ont = r] with r ~ 10 4 yr and r] = 0.4. (4) 

with t ch the characteristic pulsar age, r the characteristic luminosity decay time, E 
the spin-down luminosity, and rj the e + e~ production efficiency (our results are easily 
linearly rescaled for different values of the last parameter). For the relatively large e + e~ 
production efficiency we assume here, 77 = 0.4, the resulting total energy output in e + e~ 
E e + e - for the 10 pulsars we consider is always of the order of 10 48 erg, with the exception 
of the older pulsars J2043+2740 (for which E e+e - ~ 10 50 erg) and of J1836+5925 and 
J1057-5226 (for which E e + e - ~ 10 49 erg). Observations of pulsar wind nebulae at radio, 
X-ray and very high-energy gamma-ray energies indicate however that E e + e - > 10 48 
erg [9], with uncertainties associated to the e + e~ diffusion and magnetic fields in the 
nebulae: such large total energy outputs appear therefore observationally plausible. 

To calculate the local e + e~ flux from a source term such as the one we adopt in 
Eq. (3) we consider the following standard cosmic-ray diffusion-loss transport equation: 

9Nei ^ rl - D(E)V 2 N e - ^(b(E)N e ) = Q(E,t,f), (5) 

where D{E) = D (E/1 GeV) s is the rigidity-dependent diffusion coefficient, for 
which we assume the customary values D = 3.6 x 10 28 cm 2 /s and 5 = 0.33 
[3], and where b(E) = b E 2 is the energy loss term, which includes the dominant 
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synchrotron and inverse Compton energy loss mechanisms, with the numerical value 
b = 1.4 x 10~ 16 GeV _1 s _1 . Analytic solutions exist for the above transport equation 
with the source term we take here (see e.g. [39]), and the resulting contribution to the 
local e + e~ flux reads [3]: 



and where the maximal energy, i.e. the energy an electron or positron injected with 
arbitrarily large energy would have after a time t, is E max (t) = (bot) -1 . 

With the parameter assumptions outlined above, and with the distance estimates 
obtained in the preceding section, we obtain the results reported in Tab. 2 for the 10 
pulsars with the largest contributions to the local e + e~ flux (all other pulsars contribute 
less than 0.1% to the measured e + e _ flux at the energies corresponding to their peak 
emission). In the table we indicate, together with the pulsar name, a series of physical 
properties, including whether the object in question was discovered via the pulsar blind 
search campaign, or via radio ephemerides or whether it was an already-known gamma- 
ray pulsar (i.e. detected by EGRET), the gamma-ray spectral index, the pulsar age and 
its spin-down luminosity. The last three columns indicate the relative contribution to 
the local e + e~ flux (as measured by Fermi-LAT [2]) at an energy of E e ± = 100 GeV 
(R100) and at the energy Epeak where the flux times energy cubed peaks (Rpeak). 

We illustrate the actual e + e~ flux predictions for the ten contributing pulsars in 
Fig. 2. Red lines correspond to pre- Fermi "known" gamma-ray pulsars, including the 
Monogem (J0659+1414) and Geminga (J0633+1746) pulsars, while black lines refer to 
the new pulsars discovered with the blind search campaign on the Fermi-LAT data. 
We also show, for comparison, the Fermi-LAT measurement of the local e + e~ flux 
[2]. From the figure, it is clear that bright gamma-ray pulsars can play an important 
role as contributors to the local flux of energetic e + e~ . A potentially very significant 
contribution stems from the gamma-ray selected pulsar J202 1+4026. 

While most of the pulsars we consider here in Tab. 2 and Fig. 2 have never been 
previously studied as sources of cosmic-ray electrons and positrons, two significant and 
noteworthy exceptions (Geminga and Monogem) allow us to directly compare in detail 
with the results of recent related analyses [5-9] and [3]. In particular, our present setup 
follows exactly the same conventions as Ref. [7] and [3], so our results match with those 
presented there, modulo setting the same parameter values (most of the choices we make 
here in fact match those in [3], while a broader range of parameter choices is considered in 
[7]). Ref. [5] considers both Monogem and Geminga, although with a different diffusion 
setup (especially for the rigidity dependence 5 = 0.6 versus 5 = 0.33 adopted here) and 
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Figure 2. The spectrum of the 10 Fermi-LAT pulsars giving the largest contributions 
to the local e + e~ flux, assuming an e + e~ injection effciency r\ = 0.4 and an e + e~ 
spectral index V = 1.7 with a cutoff E cut = 1 TeV for all pulsars. Black lines refer 
to blind search gamma-ray selected pulsars, red lines to all other pulsars. The data 
points reproduce the e + e _ spectrum measured by Fermi [2] 



with a different estimate for the distance to Geminga. The range of injection spectral 
indices (1.5-2) encompasses our choice of 1.7. Our estimate for the total energy output 
for Geminga is larger than what considered in [5], although the smaller distance partly 
compensate for it. The output we estimate for Monogem is instead in the middle of 
the range envisioned in [5]. Factoring in the above mentioned differences, our results 
are in line with those presented in ref. [5] for the two considered pulsars. The spectra 
presented in ref. [6] and [9] for the Geminga pulsar appear to be slightly different from 
what we show in fig. 2. This is due to both slightly different choices for the cosmic-ray 
propagation parameters (smaller diffusion coefficient, 5 = 0.4), and to the fact that in [6] 
a time-dependence for the process of e + e~ injection in the ISM was included (leading to 
an output tail of higher energy cosmic rays), and that in [9] a fully relativistic solution 
to the diffusion-loss equation based on the Jiittner formalism is employed. In both cases, 
though, the overall predicted e + e~ output from Geminga is comparable to the present 
estimate based on gamma-ray data. Finally, Ref. [8] uses a theoretical setup which is 
very close to what we employ here, but does not single out the Geminga or Monogem 
pulsars, hence a direct comparison with our results is not directly obvious. However, 
their results on the overall e + e~ output are once again in line with what we present 
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here. 

We further elaborate on this point in Fig. 3 and 4. With the same color coding 
(black for gamma-ray selected pulsars, red for radio-selected and EGRET pulsars), we 
show in Fig. 3 the percent contribution of gamma-ray pulsars to the measured e + e~ flux 
at E e ± = 100 GeV, with the parameter space assumptions specified above, as a function 
of the pulsar age (left panel) and spin-down luminosity (right panel). The error bars 
mirror the uncertainty in the pulsar distance, as determined from our analysis based 
on gamma-ray luminosity and pulse shape. We find that at least four new Fermi-LAT 
pulsars should give a non- negligible contribution (more than 1%) at E e ± = 100 GeV. 
We thus find that the currently known gamma-ray pulsars that might contribute to the 
local e + e~ flux span a wide range of both characteristic age and spin-down luminosities, 
respectively between 10 5 and 2 x 10 6 yr and E between 0.5 x 10 34 and 10 35 erg/s. 

The relative contribution from pulsars to the e + e~ at the peak of E^ ± J (J being the 
local e + e~ flux per steradian from the given pulsar) is shown in Fig. 4, as a function of 
the pulsar distance and peak location. In the figure we indicate the horizontal errorbars 
corresponding to the uncertainty on our gamma-ray pulsar distance determination, as 
well as the resulting propagated uncertainty on the contribution at peak. As expected, 
the main pulsar e + e~ contributors lie relatively close to us, at most at a distance of 1 
kpc. The contribution from further pulsars decreases steeply with distance. 

The energy where the pulsar contribution to the e + e~ flux is expected to peak, with 
respect to the measured local e + e~ cosmic-ray flux, ranges between 100 GeV and 1 TeV. 
The location of this peak is a function of the pulsar age, that sets the maximal energy 
Emaxit) electrons and positrons can have a time t after injection. This is illustrated in 
the inset of Fig. 4, right panel, where we correlate the peak position to the pulsar age. 
The contribution from older pulsars peaks at lower frequencies, while younger pulsars 
tend to contribute at higher energies. 

4. Can Gamma-Ray Quiet Pulsars Significantly Contribute to the Local 
e + e" Flux? 

While the ratio of the gamma-ray flux from a given pulsar over the pulsar luminosity 
scales as ~ 1/-D 2 , i.e. it is inversely proportional to the square of the distance of 
the observer to the pulsar, the contribution of pulsars to the local e + e~ flux has a 
significantly more complicated dependence on distance. This stems from the fact that 
the propagation of e + e~ from the pulsar to the Earth proceeds via a diffusive process. 
Similarly, unlike the gamma-ray flux, the e + e~ flux has a critical dependence on the 
pulsar age. Here, as detailed above, we model the diffusive mechanism of charged 
particles losing energy and changing direction scattering off of the magnetic fields 
irregularities in the picture encompassed by Eq. (5). In turn, this implies the non- 
trivial and highly non-linear dependence of the pulsar contribution to the e + e~ flux on 
distance given in Eq. (6). In particular, this dependence is very sensitive to the pulsar 
age and to the assumed diffusion coefficient, via the energy-dependent diffusion radius 
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Figure 3. Percent contribution to the total e + e~ flux, as measured by Fermi-LAT 
[2], at E e + e - —100 GeV. In the left panel we show said contribution as a function of 
the pulsar age, on the right as a function of the pulsar spin-down luminosity. 



As a result of the non-trivial propagation of e + e~ in the Galaxy, it is not obvious 
which of the bright gamma-ray pulsars are the most important contributors to the 
local cosmic-ray e + e _ flux. In this section we investigate and outline the portions of 
the pulsar parameter space where pulsars are expected to contribute to the e + e~ flux, 
and the overlap of these regions with the Fermi-LAT sensitivity to the emission from 
gamma-ray pulsars. This allows us to determine which regions of the pulsar parameter 
space, important to understand the local e + e~ flux, might be within reach with further 
improvements to the Fermi-LAT sensitivity to pulsed gamma-ray emission, and which 
regions are equally important but beyond the reach of the LAT. 

In Fig. 5 we consider the parameter space plane defined by the pulsar age versus the 
pulsar spin down luminosity, for pulsars with a distance around 0.4 kpc (left) and around 
0.8 kpc (right). Pulsars within the yellow-shaded area are expected to contribute to the 
local e + e~ flux by more than 1%, for an e + e~ efficiency i] = 0.4. We also include vertical 
black lines corresponding to the pulsars with the lowest spin-down luminosity (and hence 
lowest theoretical gamma-ray flux cx L 1 (E)/D 2 ) in that given distance range. This gives 
a benchmark of the Fermi-LAT gamma-ray sensitivity for pulsars in the given distance 
range (notice that this is only to guide the eye, since the gamma-ray flux depends on 
the pulsar-dependent gamma-ray efficiency and on the geometric flux correction factor). 

We indicate, for purposes of illustration, pulsars from the ATNF catalogue (blue 
diamonds) and Fermi-LAT pulsars (black circles) with distances in the range between 
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Figure 4. Percent contribution to the total e + e~ flux, as measured by Fermi- 
LAT [2], at the peak energy for each individual pulsar. In the left panel we show 
said contribution as a function of the pulsar distance estimated from the Fermi-LAT 
measured gamma-ray flux (we also show errors on the distance estimate) , on the right 
as a function of the e + e~ flux peak energy (see also Fig. 2). The inset in the right 
panel shows the flux peak energy as a function of the pulsar age. 



0.2 and 0.6 kpc (left) and between 0.6 and 1.1 kpc (right). The corresponding yellow 
contours would of course be different for each individual pulsar at its own distance, 
but this illustrates, for instance, that the largest contribution to the local flux of e + e~ 
is expected to come from gamma-loud rather than from radio-loud pulsars (naturally, 
some gamma-ray pulsars expected to be significant contributors, such as Geminga, have 
also been detected at radio frequencies [40]). Also remarkably, regions of parameter 
space naively unaccessible to Fermi-LAT (i.e. to the left of the vertical lines) appear to 
be devoid of radio pulsars, implying that no significant e + e~ contributor is expected to 
be radio- loud and gamma-quiet. 

The upper cutoff in age for the pulsars contributing to the local e + e~ flux simply 
stems from the maximal energy cutoff as a function of energy 

B /1.4 x 10- 16 GeV-V 1 ^ /100 GeV\ 

The non-trivial dependence on pulsar age that cuts off the contribution from young 
pulsars depends, instead, on the fact that for young pulsars and for energies such that 
E -Emax = (M)~\ the diffusion radius is 

Do (W0GeV\ 5/2 

3.6 x 10 28 cm 2 s" 1 10 5 yr J \ E ) ' 



R d]S ~ 0.5 kpc 
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Figure 5. A comparison of the gamma-ray detectability of pulsars and of the predicted 
contribution to the local e + e~ flux, in the plane defined by the pulsars age versus spin 
down luminosity. In the left panel we focus on pulsars at a distance of 0.4 kpc, while 
on the right of 0.8 kpc. We also show, for reference, the location of actual pulsars, 
for distances between 0.2 and 0.6 kpc (left) and between 0.6 and 1.1 kpc (right). The 
black vertical lines correspond to the gamma-ray pulsars with the lowest spin-down 
luminosity. 



so that the contribution from a young pulsar further away than 0.5 kpc is negligible at 
100 GeV. In the picture where cosmic-ray e + e~ follow a random-walk path, this means 
that the diffusion radius of a young pulsar, at some given energy, is not large enough 
for the cosmic rays to diffuse to us. The shape of the parameter space shown in Fig. 5 
reproduces the analytic approximation outlined above. 

Fig. 6 shows the orthogonal parameter space direction of pulsar distance versus 
spin down luminosity, for pulsars with an age around 10 5 yr (left) and 10 6 yr (right). 
The black dots refer to gamma-ray pulsars, while the blue diamonds to ATNF pulsars. 
For the former, we take the distance from the present analysis, while for the latter we 
use the distances quoted in the ATNF catalogue [41]. Notice that we include all ATNF 
pulsars in the specified age ranges (0.1 — 4 x 10 5 yr in the left panel and 0.5 — 1.8 x 10 6 
yr in the right panel), irrespective of whether or not they are theoretically expected to 
produce gamma rays. In the context of the outer gap model, for instance, the condition 
for a rotation powered pulsar to produce gamma rays is that the ratio of the physical 
size of the outer gap over the ratio of the light cylinder be less than one, see e.g. Ref. [42] 
and [43]. This condition can be re-expressed in terms of the pulsar period P in seconds 
and of the pulsar magnetic field at the neutron star surface P>i 2 in units of 10 12 G as 

g = 5.5x P 26 / 21 B u 4/7 < 1. 
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Figure 6. A comparison of the gamma-ray detectability of pulsars and of the predicted 
contribution to the local e + e~ flux, in the plane defined by the pulsar distance versus 
spin down luminosity. In the left panel we consider pulsars with an age of 10 5 yr, while 
on the right of 10 6 yr. We also show, for reference, the location of actual pulsars on 
the plane, for ages in the 0.1 — 4 x 10 5 yr range (left) and in the 0.5 — 1.8 x 10 6 yr range 
(right). To guide the eye, we draw lines corresponding to the gamma-ray pulsar with 

the lowest L 1 / D 2 in the range shown, for L 1 = E (solid) and L 1 = \J 10 33 erg s" 1 x E 
(dashed). 



Not all the ATNF pulsars we show fulfill this criterion. As before, we add black lines to 
help guide the eye to the region of parameter space that should be detectable by Fermi. 
Here, we draw lines corresponding to the gamma-ray pulsar with the lowest L 7 /D 2 in 

the range shown, for L 7 — E (solid) and L 7 = J 10 33 erg s _1 x E (dashed). Objects to 
the lower-right of the lines should give (for a gamma-ray efficiency and geometric flux 
correction factor comparable to the reference pulsar) a gamma-ray flux detectable by 
Fermi-LAT. 

The shape of the parameter space where pulsars are expected to contribute to the 
local flux of e + e~ is here understood by considering the diffusion radius in the limit 
of mature enough pulsars, i.e. for energies such that E/E m3iX ~ 1. In this case, the 
diffusion radius is independent of the pulsar age, and it takes the value: 

D 1.4 x IP' 16 s \ 1/2 / 100 GeV y-* 

3.6 x 10 28 cm 2 s" 1 b ) \ E J 

The distance above sets the quadratic-exponential decay scale for the local e + e _ flux 
with the pulsar distance r, N e oc exp[— (r/i? dif j) 2 ] observed in the figure. The pulsar 



R dm ~ 2.7 kpc 
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age influences the precise location of the contour in the plane via the effect on the e + e~ 
spectrum of the maximal energy exponential cutoff. 

Fig. 6, right, shows that pulsars which have not been detected in gamma rays yet 
(those corresponding to the blue diamonds) can in principle contribute significantly to 
the cosmic-ray flux as long as (i) they have a large enough spin-down luminosity and 
(ii) they are close enough (from the figure to the right we infer a distance cutoff around 
2-3 kpc for ages on the order of a mega-year). We label some of the ATNF pulsars that 
follow in the portion of parameter space where pulsars are predicted to contribute to the 
percent level to the e + e~ flux, and that fulfill the above-mentioned gamma-ray pulsar 
outer-gap condition g < 1. We find that several of these gamma-ray quiet ATNF pulsars 
(as of the first Fermi pulsar catalogues) can potentially give significant contributions to 
the e + e~ flux. Those same objects, however, are also expected to give a potentially 
detectable gamma-ray flux, and are thus ideal targets for LAT monitoring for a pulsed 
emission. Incidentally, we notice that the age of these pulsars, t ~ 10 6 yr, indicates 
that their most significant contribution to the e + e~ flux would fall in the energy range 
around 100-200 GeV (see e.g. the inset in Fig. 4, right). 

Notice that contours of constant contribution from pulsars to the local e + e~ flux 
at 100 GeV would follow the shape of the yellow contours in Fig. 6. This means that, 
for instance, the ATNF gamma-ray quiet pulsars close to the edge and above the solid 
black line in the right panel, while potentially beyond the Fermi gamma-ray sensitivity, 
are less likely to significantly contribute to the flux of e + e~ at Earth. 

5. Summary and Conclusions 

In this paper we estimated the contribution to the local electron-positron flux from the 
gamma-ray pulsars in the first Fermi-LAT pulsar catalogue. We carried out estimates 
of the pulsar distances and their uncertainties from the Fermi gamma-ray flux and pulse 
shape. We calibrated our distance estimates with those pulsars with available alternate 
distance determinations. For a fixed choice of pulsar parameters relevant to the injection 
of electrons and positrons, we estimated the contribution of the set of known gamma- 
ray pulsars to the local high-energy cosmic-ray lepton flux. We found that ten Fermi- 
LAT pulsars can potentially contribute significantly to this flux. Of these ten potential 
contributors, six are among newly discovered pulsars, found in the first Fermi blind- 
search campaign, and have very dim or no associated radio signal. We explored the 
ranges of pulsar ages and distances where we theoretically expect gamma-ray pulsars to 
contribute most significantly to the local cosmic-ray lepton flux. We found that regions 
of parameter space exist where gamma-ray pulsars would significantly contribute to 
the e + e~ flux without producing a gamma-ray flux detectable by Fermi-LAT. However, 
known radio pulsars fall in regions of parameter space where the e + e~ contribution is 
predicted to be typically much smaller than those regions where Fermi-LAT pulsars 
exist. We also found that radio pulsars with ages around one mega-year that have 
not yet been detected in gamma rays might also contribute significantly to the local 
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e + e~ flux. In summary, we conclude that gamma-ray searches for mature and relatively 
nearby pulsars play a crucial role in understanding the origin of the local high-energy 
cosmic-ray electron and positron flux. 
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Table 1. A comparison of the distances quoted in the first Fermi-LAT pulsar catalogue 
[10] and our distance predictions based on the catalogue's gamma-ray data. 



Pulsar Name 


Catalog Distance (kpc) 


gamma-ray distance (kpc) 


J0357+32 




0.82 ± 0.16 


J 1732-31 




1.52 ± 0.31 


J1741-2054 


.38 ± .11 


1.11 ± 0.32 


J1809-2332 


1.7 ± 1 


1.45 ± 0.29 


J1836+5925 


<.8 


0.44 ± 0.4 


J2021+4026 


1.5 ± 0.45 


0.44 ± 0.06 


J0633+1746 


.25 ± .12 


0.28 ± 0.09 


J0659+1414 


.288 ± .033 


1.39 ± 0.17 


J1057-5226 


0.72 


1.02 ± 0.2 


J2043+2740 


1.8 ± .54 


3.35 ± 0.77 


J1509-5850 


2.6 ± .8 


2.87 ± 0.56 


J0631+1036 


0.75-3.62 


4.87 ± 1.25 


J1952+3252 


2 ± .5 


3.12 ± 0.44 


J0835-4510 


.287 ± .019 


0.29 ± 0.02 


J2032+4127 


1.6-3.6 


2.75 ± 0.69 


J1958+2846 




2.99 ± 0.68 


J2238+59 




3.99 ± 0.94 


J1907+06 




3.34 ± 1.03 


J1826-1256 




2.27 ± 0.48 


J1813-1246 




3.19 ± 0.69 


J 1459-60 




2.93 ± 0.76 


J1418-6058 


2.0-5.0 


3.03 ± 0.73 


J0633+0632 




2.02 ± 0.49 


J0007+7303 


1.4 ± .3 


2.41 ± 0.74 


J0534+2200 


2 ± .5 


2 ± 0.5 


J 1709-4429 


1.4-3.6 


2.5 ± 1.11 


J1048-5832 


2.71 ± .81 


2.71 ± 0.82 


J1028-5819 


2.33 ± .7 


2.33 ± 0.73 


J2021+3651 


2.1 ± 2.1 


2.1 ± 2.11 


J2229+6114 


0.8-6.5 


3.65 ± 2.85 



Pulsar Name Pulsar Type Spectral Index Age (kyr) Edot (lO' 34 erg) R100 RPeak Epcak (GeV) 



J0357+32 


BS 


1.29(0.22) 


585.00 


0.50 


8 


.74E- 


•03 


+3 
-3 


:2HE- 
.24E- 


-03 

-03 


1 


.43E- 


-02 


+4. 
-3 


06A'- 
71E- 


-03 
-03 


217l?;f 


J1732-31 


BS 


1.27(0.14) 


120.00 


13.60 


2 


.86E- 


-05 


+8 
-2 


.68E- 
.86E- 


-03 
-05 


5 


15E- 


-03 


+5. 
-7 


.04E- 
65E- 


-03 
-02 


829.84t$ : ri 


J1741-2054 


BS 


1.39(0.17) 


392.10 


0.90 


5 


.32E- 


-03 


+ 1 
-4 


14E- 
.62E- 


-02 
-03 


1 


.41E- 


-02 


+ 1 
-1 


02E- 
59E- 


-02 
-02 


303.34 


J1809-2332 


BS 


1.52(0.07) 


67.60 


43.00 


2 


.94E- 




60E- 
.94E- 


-03 
-07 


1 


.79E- 


-03 


+ 1. 
-8. 


79E- 
65E- 


-03 
-02 


1262.48±|2i;^ 


J1836+5925 


BS 


1.35(0.04) 


1800.00 


1.20 


4 


.79E- 


-02 


+ 1 
-1 


93E- 
.15E- 


-03 
-02 


4 


96E- 


-02 


+ 1. 
-2 


25E- 
07E- 


-02 
-03 


91.04±g;£S 


J2021+4026 


BS 


1.79(0.04) 


76.80 


11.60 


3 


.82E- 


-02 


+2 
-3 


89E- 
.82E- 


-02 
-02 


1 


.61E- 


-01 


+ 1. 
-1 


61E- 
.86E- 


-01 
-02 


64S.9™- 9 


J0633+1746 


EGRET 


1.08(0.02) 


342.00 


3.25 


6 


.56E- 


-02 


+(> 
-1 


15E- 
.06E- 


-03 
-02 


1 


17E- 


-01 


+ 1. 
-7. 


30E- 
01E- 


-02 
-03 


288.46i^4 


J0659+1414 


Eph 


2.37(0.5) 
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